Realizing Discontinuous Wave Functions with Renormalized Short-Range Potentials 
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We show that the most general three-parameter family of point interactions on the hne can be 
expressed as the self-adjoint local operators in terms of three Dirac's delta functions with the renor- 
malized strengths in the disappearing distances. Experimental realization of the Neumann boundary 
is discussed. 
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The discontinuity of the gradient of wave functions has 
been known since the very early days of quantum me- 
chanics. Even today, however, the admissibihty of the 
discontinuity in the wave function itself is not well recog- 
nized outside the circle of mathematical physicists. This 
is mainly due to the esoteric languages employed in for- 
mulating the subject, which tends to give an impression 
that the phenomena are far removed from the experimen- 
tally realizable settings. 

Let us focus on an example of the simplest, but non- 
trivial wave function discontinuity that can be found in 
the free particle quantum mechanics on one-dimensional 
line with a pointlikc interaction, or a point defect 
Assuming the time-reversal symmetry, the effect of the 
interaction is usually expressed in terms of the three pa- 
rameter connection condition between the wave functions 
and their derivatives at the left and the right of the loca- 
tion of the defect. Depending on the parameter value, the 
condition allows the discontinuities both of wave func- 
tion and its space-derivative. In addition, there can be 
left-right asymmetry in the connection condition at the 
defect. In the abstract mathematical treatment, how- 
ever, no specific prescription is given for the reahzation 
of this connection condition as Hamiltonian dynamics. 
Despite several attempts for the physical interpretation 
|@,||, the intuitive picture is still lacking. That is evi- 
dent, for example, in the controversy over its appellation 
of "(5' interaction" (see the discussions in |9l|ic[]). Also, 
no experiment can be conceived to check and utilize their 
effect. 

Two recent wor ks [|l]jl| on the potential models 
whose zero-range limit give the full connection condition 
go some way to address this problem. But the potentials 
obtained there are non-local and non-Hermitian except 
in the zero-range limit, thus still leave us puzzled over 
their relevance to the real-world quantum mechanics. It 
would be very useful if we can express the point defect 



in terms of a zero-range limit of well behaving functions. 

The objective of this paper is to devise just such poten- 
tial functions. It is done in terms of zero-distance limit 
of three or more Dirac's delta functions. This effectively 
gives the practical prescription to realize the wave func- 
tion discontinuity as well as its intuitive understanding, 
since, the delta function in one dimension, in tern, is re- 
alizable as a regular limit of small size potential with the 
volume integral kept constant. 

We start by defining a function using three Dirac's 
delta functions placed on the line separated by small dis- 
tances a; 



^(a:; v, u, a) — vS{x + a) + u5{x) + v5{x — a). 



(1) 



The strengths of the delta functions arc allowed to be 
varied as functions of the distance a; 



v[a), u = u{a) as a - 



0. 



(2) 



We look at the quantum-mechanical wave function '4j{x) 
under the influence of the potential ^(x; v, u, a) 



''Pix) -I- S,(x; V, u, a)'ip{x) — Eipix) 



(3) 



The solution of the Schrodinger equation in the intervals 
X G [0+,a_], where X- and x+ respectively signify the 
location infinitesimally smaller and larger than x, can be 
written in the form 



ip{x) — ijj (0) cos kx -\- iJ;'{Q^)k "'^sinfca;, 
"ip'ix) — ip'{0+) cos kx — ipiO)k sin kx, 



(4) 



where k is the wave number defined by fc = \J2E. As- 
suming that we are interested in the low energy spectra 
of Eq. (H) , we expand Eq. (^) in terms of k to take the 
leading order, and obtain 
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V'(a)=V(0) + a^'(0+), 
V''(a-) = ^'(0+). 



(5) 



With analogous treatment for x G [(— a)+, 0_], we have 

V-C-a) = V(0) - aV^'lO-), (6) 
V''((-a)+)=V'(0-). 

Note that the wave function iIj{x) itself is continuous ev- 
erywhere at this stage. Combining Eqs. (|^) and with 
the connection conditions at the location of the three 
Dirac's deltas, 

V''((-a)_) - ^'((-a)+) = -2«^(-a), (7) 
^'(0+)-V'(0-) = 2«^(0), 
^'ia+)-iP'{a^)^2vij{a), 

we obtain 

V^'(a+) - ^'((-a)_) = B{a) {^'(a) + ^(-a)} , (8) 
V-Ca) - ^(-a) = Dia) {^P'ia+) + 7A'((-a)_)} , 

where the functions B{a) and D(a) are given by 

u(a) 



B{a) = 2w(a) + 
D{a) 



1 + au{a) ' 



(9) 



2av{a) + I' 



We now take the limit a and look at the relations 
between the quantities 



■0± = lim ip(zLa), 
V-i ^ lim V''((±a)±), 

a— i-O 



(10) 



which express the possible discontinuity of the wave func- 
tion and its derivative at the location of the interaction 
which is now a single point, x = 0. 

We can choose the function v{a) and u{a) so that ei- 
ther of B{a) or D{a) becomes zero at a limit while 
keeping the other finite. One obvious choice is to keep 
v{a) and u{a) constant, that is 



v{a) = vq, u{a) — uq. 
This results in a potential function 

S{x; 2vo + uo) = lim ^(x; vo,Uo, a) 

which yields the conditions 

i^'+-Tp'_= 2(2t.o + uo)V'-, 
V'+ - = 0. 



This means 



6{x; v) = vd{x) 



(11) 



(12) 



(13) 



(14) 



which is a rather trivial result. 

Now comes the second choice, which is obtained by 
setting v{a) — > — l/(2a) to make the D{a) non-zero and 
choosing u{a) so that B{a) becomes zero; namely 



(15) 



1 


1 


2^ ^ 


2^ 
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We then have a potential 

e{x; c 
which results in 



lim £(x: , 

a^o^^ ' 2c 2a 



1 

a 



c 

2 ■ 



a) 



0, 



(16) 



(17) 



— ip- — 2c'0L 



which is what we wanted; the discontinuity of the wave 
function. 

Several remarks on the discontinuity-inducing poten- 
tial e{x; c) are in order. In our construction, the function 
is even with respect to the transposition x ^ — x. This 
fact does not seem to favour its interpretation in terms of 
"J'(a;)" as has been done in the literature. The existence 
of the third delta in the middle, which is essential to get 
the convergence, also seems to preclude the interpreta- 
tion of this function in terms of difference between two 
deltas. An important special case of £(a;; c) is obtained in 
the limit c — > oo, which results in ijj'{0) = 0. This is non 
other than the Neumann boundary condition that sepa- 
rates the system into the two regions a; > and x < 0. 
This is contrasted to the Dirichlet boundary condition, 
which is of course obtained as the v = oo limit of S{x; v). 

The self-adjoint extension theory applied to the quan- 
tum mechanics on the line with a point defect gives the 
most general connection condition for the wave function 
at the site of defect in the form |^ 



with a constraint 



-a 
-S 



aj — [35 — 1. 



(18) 



(19) 



The special choice a = 7 = — 1 gives the conditions cq. 
(|l|) (for (5 = 0) and eq. (|l|) (for (3 = Q). We therefore 
have the correspondence 
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1 

2c 



2v 
1 


1 



S{x;v), 
e{x; c). 



(20) 



In order to obtain the potential function which gives 
the full connection condition, one can work out in a sim- 
ilar fashion as before, with different choice of v{a) and 
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u{a), introducing different strengths to the deltas of left 
and right. But we take different route here, which is 
easier and, in a sense, gives us better insight into the 
meaning of the parameter values a /3 7 and S. 
The method is based on the matrix identity 



-a 
-S 



a+l 
S 
1 



1 

-S 1 



7+1 

s 

1 



(21) 



Successive application of the matrix to the vector {ip' , tp) 
can be implemented as the point interactions placed next 
to each other in disappearing distance. With the corre- 
spondence eq. ( pO| ) in mind, it is easy to convince oneself 
that the potential function 



X(x;a,/3,7,(S) 

hm ^(. + 6;^)+s(x;-) + <5(x-6; — 



(22) 



produces the desired connection condition, eq. (|l8| ) at 
X = 0. Since no singularity in the coupling is present in 
the limit 6 — > 0, one can make b arbitrarily small as long 
as it doesn't contradict with the definitions eqs.([l^) and 
(|l|). This means that one can "merge" the two deltas 
at a; = — 6 and x = h into the two deltas at a; = — a and 
X — a which constitute peripheral flank of e{x]c). One 
arrives at 



x(x;q;,/3,7,(5) 

Imi <^ 6{x + a; — — + 6{x; 

a^o 26 2a a 2a^ 



(23) 



+S{x — a; 



~26 2^' 



It is not difficult to confirm this expression with direct 
calculation similar to eqs. (0)-(§). 

The expressions eqs. ( |2ll ) - (^3|) are not valid for 6 = 0, 
in which case we resort to the expression 



-a 
-S 



1 



1 

g+l 



iMo 1 

which, in terms of the potential function, means 
x(a;;a,/3,7,(5) 

7 + 1 



(24) 



(25) 



Inn <j e{x + b; ^^^) + H^; ^) + £{x ■ 



2/3 



Combining this and eqs. (12) and (^^, and again "merg- 
ing" the neighbouring deltas where appropriate, we ob- 
tain an expression which involve five deltas; 



x(a;;a,/?,7,(5) 
5{x + 2a; 



lim 

a— >0 



f3 



+S{x; 



+S{X' 



1 



+ 



1 

f3 



2a ^ 



2 



Six 



1 

a: 

a 



7 + 1 
f^) + ^(.-2a; 



(26) 
1, 



7 + 1 



2/3a2 



20-* 



The expressions eqs. (^ 
to the special case of /3 = 
another expression 



6|) are still not applicable 
0. For this case, we use 



-a 

-7 















(27) 



1 

1 ^ 




1 



1 

_1 
p 

(1 



where p = ylaj = l/vl7l said the composite sign cor- 
responds to the case of a = T|a|. To obtain the second 
equality, eq. ( |2l| ) is applied to the both matrices in the 
RHS of the first equality. One can implement this ex- 
pression in potential form as 



X{x;a,0,j,0) 



(28) 



lim <j S{x + 2a; ±-L _ ^) + S{x + a;--±^) 
a^o I Zp 2a a Za'^ 



+5{x;-'^±- 

2 Zp a 



+5{x — a; 



-^) + 5{x~2a-~^ 



1 



a 2pa 



2 2a 



Thus, the most general connection condition around the 
point defect in one-dimensional quantum mechanics can 
be realized as singular, but renormalized zero-distance 
limit of three (or five) Dirac's delta functions. We note 
that the above expressions arc by no means unique; one 
can construct expressions with three deltas in place of 
eqs. (|6|) and (|2|) with direct method mentioned before. 
But more important than the numerical economy is the 
fact that the expressions eqs. (^2|), (|2^) and (^8|) make 
the meaning of the asymmetry for a 7^ 7 case very clear; 
one obtains different results by placing S{x; v) to the left 
or to the right of e{x; c) even in the zero-distance limit. 

We offer a numerical example to illustrate the workings 
of the realization of the wave function discontinuity. We 
use finite potential of range s in place of delta as the 
building block; 



Asix) 



2 N 

- cos i—x), 
s s 

0, 



x<\s/2\, 
X > \s/2\. 



(29) 



(30) 



Then, the potential Ea,s{x;c) defined by 

Ea,s{x]c) = (77- - ^)As{x + a) 
Zc Za 

+ (-- + ^)A,{x) + ^)As{x - a) 

a a^ 2c 2a 

becomes a good approximation to, and ultimately con- 
verge toward the e{x; c) as one takes the limit < s <C 
a ^ 0. 
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In Fig. 1, the first four eigenstates of the Schrodinger 
equation on a line x G [—L/2, L/2] with the potential eq. 
( ^ ) are shown along with the potential itself. Dirichlet 
conditions are imposed at the edge, namely, ijj{L/2) = 
^{—L/2) — 0. The value of the coupling is set to be 
c = 5. The parameters a and s, which are supposed to 
be smaller than the scale of the problem L (which we 
arbitrarily set L = 10) are chosen to be a = 0.0333i 
and s = 0.0012L. At this level of "smaU but finite" s 
and a, one can already observe the discontinuity of the 
wave function developing around x = for the second 
and the fourth states, while the first and the third states 
show "continuity" because one has V''(0) ~ for these 
even-parity states. The near degeneracy found between 
the first and the second states, and also between the third 
and the fourth states, can be thought of as the sign of the 
closeness to the Neumann limit c — oo where the system 
is divided into the two isolate identical subsystems at 
X = 0. 
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FIG. 1. (a) -(d) show the first four eigenstates of the 
particle of unit mass on the line x £ [—L/2, L/2] with the 
potential defined by eqs. (26) and (27), which is depicted in 
(e). The numbers are the energy eigenvalues. Parameters are 
c = 5, L = 10, a = 0.0333L, s = 0.0012L. The dashed lines 
show the wave functions (and the numbers in the brackets the 
eigenvalues) of the predicted a limit calculated from the 
connection condition eq. (14). (e) and (f) depict the potential 
function in different scales. 



also note that the experiments could be carried out with 
electro-magnetic and other macroscopic waves, since the 
results obtained here is applicable to larger classes of lin- 
ear wave equations other than the Schrodinger equation. 

Finally, we place our findings in broader context. 
Throughout this paper, we have kept our arguments in 
the language of elementary quantum mechanics. It is not 
difficult, however, to reformulate the problem in terms of 
Green's function In any case, we believe that our 

procedures clearly show that there are non-trivial, but 
experimentally accessible zero-range forces in one dimen- 
sional quantum mechanics other than the familiar Dirac's 
delta function. 

In hindsight, it is natural to expect the existence of 
certain singularity, or the divergence in order to obtain 
the wave function discontinuity, since the ordinary delta 
function is known to be "too week" for that purpose. For 
the theory to be still well-defined with divergent quan- 
tity, the renormalization has to be introduced. This is 
in a sense analogous to the situation in two and three 
dimensions where one encounters divergence and renor- 
malization in order to define proper point interactions 
(see, for example, and references therein). We have 
to add in rush that there is an essential dissimilarity; In 
two and three dimensions, one can define only one class 
of zero-range potential which corresponds to 6 function, 
and there is no analogue to the function e{x; c). 

It would be useful to compare the current realization of 
discontinuity-inducing potential to the earlier approach 
found in ref Q which utilizes the separable potential p^ ] 
of infinite rank, since the usual S function can be thought 
of as a separable potential of rank one. 

The function xi^', ct, (3, 7, S) belongs to a solvable class 
of quantum potential while retaining the richness of 
three parameter dependence. Already, several unusual 
physical properties are predicted for the systems with 
discontinuity- inducing potentials . Important fact to 
recall is that all imaginable potentials in one dimension 
are guaranteed to converge toward xix', a, (3, 7, 5) in zero- 
range limit. As such, it should help us unveil non-trivial 
aspects of simple quantum systems. 

In theoretical treatment of quantum system, the Neu- 
mann boundary condition appears as natural as, and at 
times, more convenient than the Dirichlet boundary. We 
hope the current work may stimulate the experimental 
design to materialize and utilize the Neumann boundary 
condition. 

We thank Prof. Izumi Tsutsui for helpful discussions. 



This example clearly shows that our procedure to re- 
alize the wave function discontinuity is not a mere math- 
ematical abstraction, but something actually realizable 
in experiments. Recent progress in the quantum device 
of nanometer scale offers a possible opportunity. We 
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